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Abstract 

We briefly review some theoretical issues involved in the direct detection of su- 
persymmetric (SUSY) dark matter. After a brief discussion of the allowed SYSY 
parameter space, we focus on the determination of the traditional neutralino detec- 
tion rates, in experiments which measure the energy of the recoiling nucleus, such as 
the coherent and spin induced rates and the dependence of the rate on the motion 
of the Earth (modulation effect). Then we examine the novel features appearing in 
directional experiments, which detect the recoiling nucleus in a given direction. Next 
we estimate the branching ratios for transitions to accessible excited nuclear levels. 
Finally we estimate the event rates leading to the atom ionization and subsequent 
detection of the outgoing electrons. 
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1 Introduction 



The combined MAXIMA-1 [1], BOOMERANG [2], DASI [3] and 
COBE/DMR Cosmic Microwave Background (CMB) observa- 
tions [4] as well as the recent WMAP data [5] imply that the 
Universe is flat [6] and that most of the matter in the Universe 
is Dark, i.e. exotic. Crudely speaking and easy to remember one 
has: 

tt b = 0.05, n CD M = 0.30, tt A = 0.65 
for the baryonic, dark matter and dark energy fractions respec- 
tively. 
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These observations, however, do not tell us anything about the 
particle nature of dark matter. This can only be accomplished 
through direct observation. Many experiments are currently un- 
der way aiming at this goal. 

Supersymmetry naturally provides candidates for the dark mat- 
ter constituents [7], [8] -[11]. In the most favored scenario of su- 
persymmetry the lightest supersymmetric particle (LSP) can be 
simply described as a Majorana fermion, a linear combination of 
the neutral components of the gauginos and higgsinos [7],[8]-[15]. 

In all calculations performed so far, the obtained event rates are 
quite low and perhaps unobservable in the near future. So one 
has to search for characteristic signatures associated with this 
reaction. Such are the modulation of the event rates with the 
motion of the Earth (modulation effect) and the correlation of 
the observed rates of directionally sensitive experiments with the 
motion of the sun [16,17]. Transitions to low energy excited nu- 
clear states have also been considered [18]. Quite recently it has 
been found that the detection of electrons, following the collision 
of the neutralino with the nucleus may offer another option [19] 
to be exploited by the experiments. 

2 The Essential Theoretical Ingredients of Direct Detection. 

The possibility of dark matter detection hinges on the nature of 
its constituents. Here we will assume that such a constituent is the 
lightest supersymmetric particle or LSP. Since this particle is ex- 
pected to be very massive, m x > 30GeV, and extremely non rel- 
ativists with average kinetic energy T « 50 keV (m x /100GeV) , 
it can be directly detected [7] -[26] mainly via the recoiling of 
a nucleus (A,Z) in elastic scattering. In this paper, however, we 
will consider alternative possibilities. 

The event rate for all such processes can be computed from the 
following ingredients: 
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(1) An effective Lagrangian at the elementary particle (quark) 
level obtained in the framework of supersymmetry as de- 
scribed , e.g., in Refs [15,23]. 

(2) A well defined procedure for transforming the amplitude ob- 
tained using the previous effective Lagrangian from the quark 
to the nucleon level, i.e. a quark model for the nucleon. This 
step is not trivial, since the obtained results depend crucially 
on the content of the nucleon in quarks other than u and d. 

(3) Nuclear matrix elements [27] — [29], [18], [19], obtained with 
as reliable as possible many body nuclear wave functions. 
Fortunately in the most studied case of the scalar coupling 
the situation is quite simple, since then one needs only the 
nuclear form factor. Some progress has also been made in ob- 
taining reliable static spin matrix elements and spin response 
functions [29], [18]. 

The calculation of this cross section has become pretty standard. 
One starts with representative input in the restricted SUSY pa- 
rameter space as described in the literature for the scalar inter- 
action [13,15] (see also Arnowitt and Dutta [24]). 

It is worth exploiting the contribution of the axial current in the 
direct neutralino detection, since, in addition, it may populate 
excited nuclear states, if they happen to be so low in energy that 
they become accessible to the low energy neutralinos [18]. Models 
which can lead to detectable spin cross sections have recently 
been proposed [20] [21] [22]. 

Once the LSP-nucleon cross section is known, the LSP-nucleus 
cross section can be obtained. The differential cross section with 
respect to the energy transfer Q for a given LSP velocity v can 
be cast in the form 

du — — 
da{u, v) = bv y ^sF\u) + S spm F n (w)] (1) 

where we have used a dimensionless variable u, proportional to Q, 
which is found convenient for handling the nuclear form factors 
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Table 1 

The static spin matrix elements for various nuclei. For light nuclei the calculations 
are from Divari et al (see text) . For 127 1 the results are from Ressel and Dean (see 
text) (*) and the Jyvaskyla-Ioannina collaboration (private communication) (**). 
For 209 Pb they were obtained previously (see text). 





19 F 


29 Si 


23 Na 


127 j* 


127 j** 


207p 5 + 


[tto(0)] 2 


2.610 


0.207 


0.477 


3.293 


1.488 


0.305 


[^i(O)] 2 


2.807 


0.219 


0.346 


1.220 


1.513 


0.231 


fio(o)JMo) 


2.707 


-0.213 


0.406 


2.008 


1.501 


-0.266 


IHh 


2.91 


-0.50 


2.22 








l^exp 


2.62 


-0.56 


2.22 








H th (spin) 

flexp 


0.91 


0.99 


0.57 









[26] F(u) , F n (u), namely u = ^ , Q « 40 x MeV. fi r is 
the reduced LSP-nucleus mass and b is (the harmonic oscillator) 
nuclear size parameter. Furthermore 



2 spin 



Z spill, > > 



spin 



f°A\2 



3(1 + ft) 



S(u)(2) 



cr^Flo and a* ^ are the nucleon cross-sections associated with the 

PiX. fiX 

spin and the scalar interactions respectively and 

5(1*) = [(^o(O)) 2 ^ + 2^0(0)^(0)^4 + ^i(0)) 2 ](3) 



The definition of the spin response functions F^, with z, j = 0, 1 
isospin indices, can be found elsewhere [29]. 

Some static spin matrix elements [29], [27], [26] for some nuclei 
of interest are given in table 1 
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3 Rates 



The differential non directional rate can be written as: 

dRundir = — t~~ da(u, v)\v\ (4) 



Where p(0) = 0.3GeVy cm 3 is the LSP density in our vicinity, m 
is the detector mass, m x is the LSP mass and da(u, v) was given 
above. 

The directional differential rate, in the direction e of the recoiling 
nucleus, is given by : 



dRdir = ^ ^ —. M"0.e @(v.e) — da(u,v) (5) 

m x ArriN ztt ' 

\fu 1 , x > 

fi r vbV2 0,x<0 



The LSP is characterized by a velocity distribution. For a given 
velocity distribution f(^'), with respect to the center of the galaxy, 
One can find the velocity distribution in the lab frame f(v, ve) 
by writing v = v+ ve , ve=Vq+ V\. vq is the sun's 
velocity (around the center of the galaxy), which coincides with 
the parameter of the Maxwellian distribution, and V\ the Earth's 
velocity (around the sun). Thus the above expressions for the 
rates must be folded with the LSP velocity distribution. We will 
distinguish two possibilities: 

(1) The direction of the recoiling nucleus is not observed. 
The non-directional differential rate is now given by: 

,dR un dir\ ,dR x p(0) m r—^-dH 
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where 



,gE. r \v\ da(u,v) 7 o 

W = /7P) / ^^ ) T d " (7) 

(2) The direction e of the recoiling nucleus is observed. 

In this case the directional differential rate is given as above 
with: 

d£ f v.e ejv.e) da(u,v) 

— m v.e)d 6 v 

To obtain the total rates one must integrate the two previous 
expressions over the energy transfer from Q m i n determined by 
the detector energy cutoff to Q max determined by the maximum 
LSP velocity (escape velocity, put in by hand in the Maxwellian 
distribution), i.e. v esc = 2.84 vq, vq = 229 Km/ s. 



4 Results 



We will specialize the above results in the following cases: 



4-1 Non directional unmodulated rates 



Ignoring the motion of the Earth the total non directional rate is 
given by 

R = Rt{a,Q mm ) (9) 



m x o Am p Hr{P) p,x 
where t is the modification of the total rate due to the folding 
and nuclear structure effects, t depends on Q m i n , i.e. the energy 
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Rbar Rbar 




Fig. 1. On top:The quantity R, « event rate for Q m in = 0, associated with the spin 
contribution in the case of the A = 19 system (for the definition of the parameters 
see text), bottom: The event rate, associated with the spin contribution in the case 
of the A = 127 system (for notation see our earlier work [18]). 

transfer cutoff imposed by the detector and the parameter a intro- 
duced above. All SUSY parameters, except the LSP mass, have 
been absorbed in R. 

Via Eq. (9) we can, if we wish, extract the nucleon cross section 
from the data. For most of the allowed parameter space the ob- 
tained results are undetectable. As it has already been mentioned 
it is possible to obtain detectable rates in the case of the coherent 
mode. Such results have, e.g. been obtained by Cerdeno et al [30] 
with non universal set of parameters and the Florida group [31]. 

For the target 19 F are shown in Fig. 1 (top), while for 127 1 the 
corresponding results are shown in Fig. 1. 



4-2 Modulated Rates 



If the effects of the motion of the Earth around the sun are in- 
cluded, the total non directional rate is given by 

R = Rt[(l + h{a, Q min )cosa)} (10) 



7 



Table 2 

The parameters t, h, re, h m and a m for the isotropic Gaussian velocity distribution 
and Qmin = 0. The results presented are associated with the spin contribution, but 
those for the coherent mode are similar. The results shown are for the light systems. 
For intermediate and heavy nuclei there is a dependence on the LSP mass. +x is 
radially out of the galaxy (0 = it/2, $ = 0), +z is in the sun's direction of motion 
(6 = 0) and +y is vertical to the plane of the galaxy (0 = tt/2, $ = ir/2) so that 
(x, y, x) is right-handed. a m = 0, 1/2, 1, 3/2 means that the maximum occurs on the 
2nd of June, September, December and March respectively. 



type t h dir re h m a Ti 



+z 0.0068 0.227 1 

dir +(-)x 0.080 0.272 3/2(1) 

+(-)y 0.080 0.210 (1) 

-z 0.395 0.060 



all 1.00 

all 0.02 



with h the modulation amplitude and a is the phase of the Earth, 
which is zero around June 2nd. The modulation amplitude would 
be an excellent signal in discriminating against background, but 
unfortunately it is very small, less than two per cent (see table 
2). Furthermore for intermediate and heavy nuclei, it can even 
change sign for sufficiently heavy LSP. So in our opinion a better 
signature is provided by directional experiments, which measure 
the direction of the recoiling nucleus. 

4-3 Directional Rates. 

Since the sun is moving around the galaxy in a directional exper- 
iment, i.e. one in which the direction of the recoiling nucleus is 
observed, one expects a strong correlation of the event rate with 
the motion of the sun. The directional rate can be written as: 



Rdir = -^ L R[^ + h m cos(a - a m tt)] (11) 

Z7T 



8 



= — R t[l + h m cos(a — a m tt)} 

where h m is the modulation, a m is the shift in the phase of the 
Earth a and k/(2tt) is the reduction factor of the unmodulated 
directional rate relative to the non-directional one. The param- 
eters k , h m , a m depend on the direction of observation: The 
above parameters are shown in Table 2 parameter t^ r for a typ- 
ical LSP mass 100 GeV is shown in for the targets A = 19. For 
heavier targets the depend on the LSP mass [16]. 

4-4 Rates to excited states 

Transitions to excited states are possible only for nuclear sys- 
tems characterized by excited states at sufficiently low energies 
with quantum numbers, which allow for Gamow- Teller transi- 
tions. One such system is 127 1, which, fortunately, can serve as a 
target for the recoil experiment. 

This nucleus has a ground state 5/2 + and a first excited state a 
7/2 + at 57.6keV. As it has already been mentioned it is a popular 
target for dark matter detection. As a result the structure of its 
ground state has been studied theoretically by a lot of groups 
(for references see [18]). We find tt 2 = = fio^i = 0.164, 0.312 
for the ground state and the excited state respectively. 

In presenting our results it is advantageous to compute the branch- 
ing ratio. In addition to factoring out most of the uncertainties 
connected with the SUSY parameters and the structure of the 
nucleon, we expect the ratio of the two spin matrix elements to 
be more reliable than their absolute values. Taking the ratio of 
the static spin matrix elements to be 1.90 and assuming that the 
spin response functions are identical, we calculated the branch- 
ing ratio , which is exhibited in Figs 2 and 3. We notice that the 
dependence on Q mm is quite mild. From Figs 2 and 3 we notice 
that the relative modulation is higher when the phase space is 
restricted by Q m i n and E exc at the expense, of course, of the total 
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Fig. 2. On the left we show he ratio of the rate to the excited state divided by 
that of the ground state for 127 1 assuming that the static spin matrix element of 
the transition from the ground to the excited state is a factor of 1.902 larger than 
that involving the ground state, but the spin response functions are the same. Next 
to it we show the modulation amplitudes for the ground and the excited states 
respectively. The results were obtained for no threshold cut off (Q m in = 0). 




Fig. 3. The the same as in Fig. 2 for a lower energy cutoff of Q m in = WkeV. 

number of counts. 



4-5 Detection of ionization electrons 



The differential cross section for the LSP nucleus scattering lead- 
ing to emission of electrons in the case of non relativistic neu- 
tralinos takes the form: 



T. Vnt [fa m (k)] 2 2n6(T x + e nf - T - ^- - (P °" k " q)2 ) 

where v is the oncoming LSP velocity, M is the invariant ampli- 
tude, known from the standard neutralino nucleus cross section, 
T and k are the kinetic energy and momentum of the outgoing 
electron, q is the momentum transferred to the nucleus and e n e is 
the binding energy of the initial electron. (f) n £ m (k) is the fourier 
transform of the bound electron wave function, i.e its wave func- 
tion in momentum space. p n e is the probability of finding the 
electron in the n£ orbit [19]. 
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Fig. 4. The cross section for ionization divided by the standard LSP nucleus cross 
section function of the threshold energy. The full line results by including the 

Is, 2s and 2p orbitals, while the dashed line is obtained considering the Is orbital 
alone 

In order to avoid any complications arising from questions re- 
garding the allowed SUSY parameter space, we will present our 
results normalized to the cross section of the standard neutralino 
nucleus cross section, which is fairly well known. One then can 
perform the needed integrals to obtain the total cross sections. 

In addition, of course, one must convolute the above expression 
with the velocity distribution to obtain both the differential cross 
section as well the total cross section as a function of the neu- 
tralino mass. We will not, however, address these issues here. 

We will now apply the above formalism for a typical case, namely 
m x =100 GeV and T x =< T x >= 50 keV, which is the average 
energy for this mass. We will also use as a target 19 F, which is 
considered in the standard nuclear recoil experiments. The ob- 
tained results are essentially identical to those for 20 Ne, which 
is a popular gaseous TPC counter currently being considered for 
detection of low energy neutrinos produced in triton decay. 

The binding energies employed [32] are 

e 0s = -0.870 , e 2s = -0.048 , e 2p = -0.021 

Our results were obtained considering one electron per atom and 
weighing each orbit with the probabilities^ = (2/10,2/10,6/10) 
in the above order. The total cross section for this process divided 
by that of the standard LSP-nucleus cross section is shown in Fig. 
4, without the inclusion of the nuclear form factor on the left and 
using appropriate form factor [29] on the right. From this plot we 
see that the introduction of the nuclear form factor decreases the 



branching ratio by a factor of about 4. This reduction maybe 
somewhat less, if the convolution with the velocity distribution 
is included. 

5 Conclusions 

In the present paper we have discussed the parameters, which de- 
scribe the event rates for direct detection of SUSY dark matter. 
In the coherent case, only in a small segment of the allowed pa- 
rameter space the rates are above the present experimental goals 
[13,15,24], which, of course, may be improved by two or three 
orders of magnitude in the planned experiments [33]- [37]. In the 
case of the spin contribution only in models with large higgsino 
components of the LSP one can obtain rates, which may be de- 
tectable, but in this case, except in special models, the bound 
on the relic LSP abundance may be violated. Thus in both cases 
the expected rates are small. Thus one feels compelled to look 
for characteristic experimental signatures for background reduc- 
tion, such as correlation of the event rates with the motion of 
the Earth (modulation effect) and angular correlation of the di- 
rectional rates with the direction of motion of the sun on top of 
their seasonal modulation. Such experiments are currently under 
way, like the UKDMC DRIFT PROJECT experiment [38], the 
Micro-TPC Detector of the Kyoto-Tokyo collaboration [39] and 
the TOKYO experiment [40]. 

The relative parameters t and h in the case for light nuclear tar- 
gets are essentially independent of the LSP mass, but they depend 
on the energy cutoff, Q m in- For Q m i n = they are exhibited in 
Table 2. They are essentially the same for both the coherent and 
the spin modes. For intermediate and heavy nuclei they depend 
on the LSP mass [16]. 

In the case of the directional rates it is instructive to first sum- 
marize our results regarding the reduction of the directional rate 
compared to the usual rate, given by k/(2tt). The factors k de- 
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pend, of course, on the angles of observation (see Table 2). Second 
we should emphasize the importance of the modulation of the di- 
rectional rates. In the favored direction the modulation is not 
very large, but still it is three times larger compared to that of 
the non directional case. In the plane perpendicular to the sun's 
motion the modulation is quite large (see Table 2). 

Coming to transitions to excited states we believe that branching 
ratios of the size obtained here for 127 1 are very encouraging to the 
experiments aiming at 7 ray detection, following the de-excitation 
of the nucleus. 

Regarding the detection of the emitted electrons in the LSP- 
nucleus collision we find that, even though the distribution peaks 
at low energies, there remains substantial strength above 0.2 keV, 
which is the threshold energy of a Micromegas detector, like the 
one recently [41] proposed. We should emphasize that this region 
is below the threshold of the nuclear recoils. 

We thus hope that, in spite of the experimental difficulties, some 
of the above signatures can be exploited by the experimentalists. 
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